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Abstract 

There are at least two different putative steady state solutions for current across noncollinear 
magnetic multilayers; one has a discontinuity in the spin current at the interface the other is 
continuous. We compare the resistance of the two and find the solution with the continuous spin 
currents is lower. By using the entropic principle we can state that this solution is a better estimate 
of the resistance for a noncollinear magnetic multilayer 
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When current is driven across a multilayered structure in which the magnetic layers 
are not collinear one has to consider how the spin polarized current is transmitted between 
layers. In most treatments one considers the components of the spin current transverse to the 
magnetization of a layer is either absorbed or reflected at the interface upon entering the layer 
and that this is the origin of spin torque created by spin currents in magnetic multilayers 1 . 
This scenario does not envisage the existence of components of the spin accumulation that 
are transverse to the magnetization of the layer into which they are injected by the charge 
current, but it does assume a steady state of the current is reached, e.g., one posits that 
there is no time dependence of the current densities. An alternate proposal is to inject 
transverse spin currents through interfaces that propagate in transition-metal ferromagnets, 
e.g., in Co up to 3 nm from the interface, and to have the spatial variation of the spin 
currents balanced by transverse spin accumulation thereby achieving steady state.— In this 
case the spin currents create an accumulation which acts as a torque on the background 
rather than directly suppling angular momentum to the background magnetization.— When 
the current drives a system out of equilibrium it is not immediately obvious which of the 
two scattering processes around the interfaces better describe the steady state transport 
across a multilayered structure. To resolve this question we use the variational principle of 
transport that says that: "in the steady state, the currents are such that the production 
of entropy has the largest value consistent with its subsequent conservation."— We find in 
the steady state transverse spin currents with their attendant accumulations do exist in 
magnetic layers in the immediate region about the interface, i.e., within several nm. Their 
existence determines the amount of spin torque created on a magnetic layer for a fixed 
current; therefore it changes one's predictions of the critical current needed to switch the 
magnetization in a multilayered structure.— 

The experiments that have demonstrated current induced magnetization reversal have 
been done on nanopillars that consist of two magnetic layers separated by a nonmagnetic 
spacer to magnetically isolate the layers, but otherwise transmit the spin polarized current 
from one layer to the other, and nonmagnetic leads£ To resolve the issue of how the spin 
polarized current created in one layer is transmitted to another noncollinear layer we solve 
for the transport across two magnetic layers and assume the nonmagnetic spacer is thin 
enough to be neglected, i.e., when one drives a spin polarized current across the trilayer one 
uses spacer layer thicknesses small compared to the spin diffusion length so that one can 
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coherently transmit spin information from one ferromagnetic layer to the other;— the analysis 
can be repeated with a finite spacer without changing our conclusions. We solve the Boltz- 
mann equation for current across the layers by using the interface scattering matrix to relate 
the distribution functions across their common interface to self-consistently determine the 
accumulations and currents in each layer, e.g., see Waintal and Bauer.— The additional resis- 
tivity due to the noncollinearity of the layers is found from the change in chemical potential 
that develops across the structure, e.g., see Valet and Fert.— While this is proportional to 
the angle between the magnetization of the layers, it also depends on whether the scattering 
at the interface is able to create spin flips; recently we have shown that they are induced 
by currents in noncollinear magnetic multilayers.— By determining the change in chemical 
potential across the two layered structure we use the variational principle, that holds that 
the electrical resistivity is a minimum in the steady state,— to discern which scattering at 
the interface achieves the lowest resistance and therefore is the best description for a steady 
state spin current in noncollinear multilayers. 

We consider two bands crossing the Fermi surface of a ferro magnet, a majority spin band 
and one minority; this is appropriate for example for Co where the two additional minority 
bands could be included but do not add much to the transport across the interface when 
the layers are noncollinear. Our Boltzmann equations for two band are,— 

Up ~ Up)) ~ — ■ Up ~ U P ')) , W 

~ ■ (f P ± - <£» , (2) 
p 

where p,p' = M for majority or m for minority, J is the exchange splitting between the 
state on the Fermi surface and its partner with the same momentum and opposite spin,— 
and r p and r s f represent spin dependent but non-spin-flip and spin-flip relaxation times. 
Also, Tp represent the transverse relaxation time; in our calculation we take t' p = r p as it 
does not affect the point we are making in this paper. The distribution functions f p (k, x) 
and fp{k,x) are functions of momentum and position; the angular brackets on the right 
hand side of Eqs.(l) and (2) denote averages over the Fermi surface. In our notation each 
band crossing the Fermi surface is represented by a spinor distribution function; these are 
statistical density matrices which characterize the out of equilibrium state of the current. 



df P 



V ~ p 'Jx~~ eE & v P 



de 



■ j p 

dx 



± 



3 



The first equation refers to the longitudinal or spin diagonal component on the distribution 
function on the Fermi surface, and the current induced coherences with states off the Fermi 
surface are given by the second equation. We do not write an equation for the excited state 
itself as it is not occupied in equilibrium; this term corresponds to the second longitudinal 
or diagonal component of the spinor. As we limit our current analysis to linear response in 
the electric field we do not consider it further here. 

As we consider two bands, p = M, m, there is a spinor for each sheet of the Fermi surface. 
For specificity we rewrite the distribution functions as, 

I'm = I'mm, /m = I'mm', I'm = fM'M, and f m = f mm , = f m ' m , f m = fmm'i (3) 

where the prime denotes an excited state. Note that although there is scattering on the 
Fermi surface to produce finite values for /a/™, fmM we do not consider these components 
because their contribution to the current cannot be calculated in the Boltzmann approach 
as there is no unique velocity associated with them.^ As we have coherences with excited 
states of opposite spin, and not only with states of opposite spin on the Fermi surface, we 
have in linear response six components of the distribution functions to keep track of in each 
layer which is nearly twice the number when one limits oneself to scattering on the Fermi 
surface. The boundary condition for these functions at the interface are — 

fqq'(® + ) = /jTpp'qq'fm'i^) + R PP'll' fpp' ( 0+ ) ' ( 4 ) 

pp' pp' 

= Ev«'/; P '(o + ) + J2 R pp'<*'fpA°^> ( 5 ) 

pp' pp' 

where pp', qq' = MM, MM' , M'M, m'm, mm' or mm. The transmission and reflection coef- 
ficients ( probabilities) are 

Tpp'=^qq' Tppiqqi tpq * tpiqi , 

Rpp'^qqi — Rpp'qq' = T pq * T plql , (6) 

where p,p',q,q' = M, M' ,m' ,m , and the amplitudes tp q ,r pq should be determined from 
ab-initio calculations.— These coefficients have been determined for interfaces between fer- 
romagnetic and normal (F/N) metals for states that exist on the Fermi surface when the 
system is in equilibrium; in this case states of opposite spin are coherent in the normal 
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layer but not in the ferromagnet when one uses a single particle description for the electron 
states. For this reason distribution functions in the ferromagnetic layer f qq i only with q = q' 
exist when the components are referred to the axis of the magnetization of the layer. We 
have recently shown that to calculate the conductivity of noncollinear magnetic multilayers 
by using a local layer-by-layer approach it is necessary to include additional scattering at 
the F/N interfaces; this produces so called current induced coherences in the ferromagnetic 
layers, i.e., components of the distribution function f qq > with q ^ g'.— In contrast to the 
case when one neglects these current induced coherences we have 16 components of to 
consider. 

To discuss transport across a noncollinear F/N/F trilayer we write the transmission and 
reflection amplitudes and coefficients in terms of the F/N scattering coefficients discussed 
above. We will make the simplification of neglecting the additional scattering at the F/N 
interfaces arising from differences in the band structure between the ferromagnetic and 
normal layers, i.e., we overlook differences between k n and kM/m, but we do consider fc M 7^ 
oversight in no way alters the primary conclusion we arrive at. When the thickness 
of the normal spacer is small compared to the spin diffusion length there is no loss of spin 
information and one can write the scattering amplitudes and coefficients connecting the 
ferromagnetic layers in terms of those for the F/N interfaces.— We use the variational principle 
cited above to determine whether non spin-flip or spin-flip scattering at the interface yields 
the best estimate for the resistivity. To calculate the angular dependence of the reflection 
and transmission due to this perturbation we use the spin rotation matrices to refer the 
states quantized along the magnetization of one F layer in terms of states quantized along 
the magnetization of the other F layer. In the limit of neglecting the non-spin-flip scattering 
at the individual F/N interfaces in which the electron stays on the same sheet of the Fermi 
surface we find for the majority state the transmission amplitudes are 



tMM = tM'M' = COs(0/2), 



(7) 



turn 1 = tu'M = iAsm9sin(8/2), 



(8) 



thlm' — tht'm — 0, 



(9) 
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2&m(1 — A sin 9) . . ,„,„,, 

iMm = *MW = ^ , , . > hl(fl/2)), (10) 

where &m, &m are the Fermi momenta of the majority and minority bands, 9 the angle be- 
tween the magnetic layers, and A parameterizes the amplitude arising from current induced 
spin-flip scattering at the F/N interfaces; A = corresponds to no spin flip scattering which 
is the conventional assumption. The half angle terms arise from the spin rotation matrices 
while the sin 9 comes from the magnitude of the current induced spin-flip scattering.— For the 
transmission amplitudes of the minority state we replace M — > m, m — > M, ku — > k m , and 
A —>■ B . When we include the scattering at the F/N interfaces that comes from ku/m 7^ k n , 
thiM and t mm are less than 1 when 9 = 0, but this does not change the conclusions we arrive 
at here.— 

For the reflection amplitudes we write for the majority state 

r M M = r M 'M> = t — — r^ 1 - Asin0)(l - cos0), (11) 

% + Km 

. ku — k m . 2 Q /1 oN 

r M M' = r M 'M = -It — -—7 A sin 9, (12) 

^Mm 1 = TMm = fM'm' = r M'm = 0, (13) 

while those for the minority state are obtained by the following replacements M — >• m,m — > 
M, ku k m , and A —> B .The unitary of scattering matrix can be confirmed as Yl qq ' Tpp> qq i + 
J2 qq ' Rpp'qq' — 1- Based on a specific mechanism, e.g., the current induced spin-flip suggested 
in Ref.(8), the values for A and B can be determined, but their precise values are unimportant 
to our argument; only their existence is crucial. 

To determine whether current induced spin-flip scattering is needed to achieve steady 
state current conditions we compared the resistance for A = B = to the resistance one finds 
for A = B ^ ; as our result is not sensitive to the values for them, in our calculations we 
consider two set of values A = B = 0.05 and 0.5. Also we used the Fermi momenta ku-, k m for 
Cc 12 so that ^~^ m — 0.4. We solved Eqs.(l) and (2) for the out of equilibrium spin transport 
across two noncollinear magnetic Co layers by using the above transmission and reflection 
amplitudes and the coefficients Eqs. (jBJ) to match the out of equilibrium distribution functions 
across the F/N/F interface. In Figs. 1 through 3 we compare the spin accumulations 
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FIG. 1: The spin accumulations for two Co layers for 6 = 90° calculated by solving Eqs.(l) and 
(2) and using the transmission and reflection amplitudes and coefficents given in the text. In 
(a) we consider no spin-flip scattering the interface, A = B = 0; while there is longitudinal spin 
accumulation m z (parallel to the local magnetization) the transverse (to the local magnetization) is 
zero, i.e., there is no injection of transverse spin currents in the absence of the interfacial spin-flip 
scattering. In (b) we have a very small amount of spin-flip scattering A = B = 0.05 and find 
a small amount of transverse spin accumulation m x and m y while the longitudinal is practically 
unchanged. For spin-flip scattering ten times stronger A = B = 0.5 the transverse accumulation is 
proportionately larger, but the range is unchanged (see Fig. 2),; the longitudinal remains unchanged. 

and currents across the interfacial region for the three sets of A = B , for 9 = 90°. For 
A = B = there is no transverse accumulation nor current (j + = j x + ij y ); only the 
longitudinal component of the spin current j z exists in the ferromagnetic layers as expected. 
For A = B 7^ the transverse spin current exists and is independent of the amount of 
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FIG. 2: The modulus of the transverse accumulation m + = m x + im y for the two values of the 
spin- flip scattering considered in Fig.l. Note the scale change for the accumulation by a factor of 
ten; otherwise the two accumulations are very similar. This underscores that the length scale of 
the transverse spin accumulation is determined by the exchange parameter J in the equation of 
motion for the transverse distribution function (see Eq.(2)), and not by the amount of spin- flip 
scattering at the interfaces. 

spin-flip scattering at the interface; only the amount of transverse accumulation depends 
on the A and B. From Fig. 2 we see that the length scale of the transverse accumulation 
(m + = m x + irtiy) does not depend on A = B; rather it is controlled by the exchange 
splitting J in Eq.(2) for the transverse distribution function. Note that the vertical scales 
for the accumulation in Figs. 2 differ by a factor of ten, so that the transverse accumulation 
is nearly linear dependent on the spin-flip scattering parameter A. 

The additional resistance due to noncollinearity of layers is found from the additional 
voltage drop AV of the electrochemical potential; see Valet and Fert.^ In the present notation 
this electrochemical potential \i e is denned as (/mm) + (fmm) where the brackets denote 
averages over the Fermi surface, and the additional resistance is 

**-? -if (§£-*)*■ (14) 

Je Je J -oo \ e udj J 

where j e = ^^-e 2 E X (v 2 M + v^), A is the mean free path in the bulk of the metal, and for 
Co the ratio of the Fermi velocities is Vm : v m = 1 : 0.78 . While the integral is over the 
entire structure we can cut it off when the chemical potential is independent of position. 
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FIG. 3: The spin currents across two noncollinear Co layers with = 90 calculated from the 
solutions for the distribution functions(see Eqs.(l) and (2)); j 2 is the component parallel to the 
local magnetization, and j+ = j x + ij y is transverse. The longitudinal and transverse spin currents 
are unchanged as we vary the amount of spin-flip scattering at the interface, albeit for A = B = 
there is no injection of transverse spin current. 

As the spin diffusion length for a typical 3d transition- metal, e.g., Co, is of the order of 60 
nm we have limited the integration to 90 nm about the interface. We have calculated the 
resistances for the three sets of A and B as a function of angle and present our results in 
Fig. 4. While the lowering of the resistance is proportional to the amount of current induced 
spin-flip scattering (A and B) the signal conclusion is that as long as A and B are nonzero 
the resistance is always lower than for A = B = 0. We conclude that when the spin current 
is discontinuous at the F/N interfaces, i.e., A = B = the resistance is always larger than 
when the current is continuous (A = B ^ 0). By invoking the principle that the electrical 
resistance is a minimum in the steady state of the current, we conclude that the spin-flip 
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FIG. 4: The resistance across two Co layers, each 90 nm, as determined from Eq.(|14j) as a function 
of the angle between their magnetizations. As we neglect the non-spin- flip scattering at the F/N 
interfaces the resistance at 9 = is zero and that for 9 = 180° is independent of the spin-flip 
scattering A = B at the interface; all resistances are normalized to that at 9 = 180°. The salient 
feature is that the presence of spin-flip scattering at the interface invariably lowers the resistance. 

interface scattering is necessary to achieve steady state conditions. 
The resistance shown in Fig0] can be written as 

R = Ra + Re*f(0), (15) 

where Rq is the resistance coming from the band mismatch at the F/N interfaces, i.e., 
k>M/m 7^ k n ; we have set this to zero in our calculation. For A = B = we are able to fit our 
results to the form suggested by^ 

1- cos 2 0/2 

/W -l + X cos^/2' (16) 
and we find \ = 3-4. For A = B = 0.05, x = 3.7; while for A = B = 0.5 this form does not 
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provide a suitable fit to the angular dependence of the resistance. However when we use 



1- cos 0/2 
/W -l + X cos£/2' 

we are able to obtain an excellent fit with \ = 1-07 for A — B — 0.5. 

In equilibrium only states up to the Fermi level are occupied and the scattering at the 

F/N interfaces is characterized by A = B = 0. When the current is switched on the current 

induced spin-flip scattering at the interfaces induces coherences between states on the Fermi 

surface and excited states of opposite spin. This injection of transverse components of the 

spin current is offset by the interference coming from the individual k states on the Fermi 

surface,— so that after a characteristic length X tr the current is parallel to the background 

magnetization; from Fig.3 we can see that Xt r — 3nm£ While in linear response we are 

only interested in the coherences f pp r between the spin states on the Fermi surface p and 

the excited states of opposite spin p', it follows from the inequality for density matrices 

representing statistical mixtures of states that 



fppfp'p' > \fpp'\ , (18) 

i.e., the presence of coherences indicate that there is an occupancy of the excited state 
when the system is out of equilibrium; albeit it is a higher order effect. As the transverse 
components of the spin current disappear close to the interface (within 3-5 nrn^) the time 
required to achieve this "equilibrium" is short compared to the spin-flip time that equilibrates 
the longitudinal components of the spin current. 

The existence of finite current induced coherences f vv > close to the interface give rise to 
the transverse spin accumulations found in this and previous studies; it should be stressed 
that our self consistent solutions of the Boltzmann equation do not mandate transverse spin 
accumulation, rather they are found to be necessary concomitant to achieve a steady state 
in noncollinear magnetic structures.— It is in this sense that our calculations indicate that 
spin-flip scattering at the interfaces is necessary to achieve steady state current conditions in 
a noncollinear multilayered structure. Finally, the time necessary for creating the transverse 
spin accumulation is short compared to that for the longitudinal; this is controlled by the 
spin-flip time which for Co is of the order of 10 -12 sec. Both times are short compared 
to the time it takes for the background magnetization to move, ~ 10~ 9 sec; therefore one 
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can achieve steady state currents in noncollinear magnetic structures prior to magnetization 
reversal provided one allows for transverse spin accumulation in the regions about the in- 
terfaces. Conversely, if one does not allow this one cannot achieve steady state conditions 
in noncollinear magnetic structures. 

In conclusion, in the calculations whose results we presented we have taken into account 
differences in the band structures in the ferromagnetic layers, i.e., ku ^ k m . In the con- 
ventional approach, A = B = 0, there is a discontinuity in the spin current at the N/F 
interfaces, while in our approach, A = B ^ 0, there is none. Therefore the origin of the 
discontinuity cannot be the band mismatch per se; rather the key difference between the 
two scenarios is that we account for the coherence between states of opposite spin in the 
ferromagnetic layers that is lost in effective single electron treatments. When we account 
for these coherences in our model with spin split bands in the ferromagnetic layers the spin 
currents are close to those found in previous calculations based on unsplit free electron bands 
(which maintain coherences throughout the multilayer). 

We would like to acknowledge very helpful discussions with Jerry Percus on the use of the 
resistivity minimum principle for steady state currents, and to Shufeng Zhang on general 
transport theory. This work done under a grant from the National Science Foundation(Grant 
DMR 0131883). 
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